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Numerical  simulations  were  performed  in  support  of  the  SRI  small-scale 
laboratory  experiments  of  deep-based  structures.  The  experiments,  in  which 
small  cylinders  of  rock  simulant  SRI  RMG  2C2  witli  lined  tunnels  were  loaded 
statically  and  dynamically  to  represent  deep-based  structure  response,  were 
investigated  analytically  by  use  of  the  CRT/NONSAP  nonlinear  finite  element 
computer  code.  The  computer  code,  which  employed  refined  material  models  for 
the  rock  simulant  was  used  to  simulate  the  SRI  test  procedure  for  both  static 


DD  1473  E'l.TiON  O'  I  NOv  AT  ilOUiHFil 


UNCLASSIFIED 

ST  C.  u  III  T  *  C  •_  A  t  1  A  '  ■  IN  O  »  T  •<  S  >-  A  ,1  l>l  ..  1.  •  r~ 


! 


UNCLASSIFIED 


*ecu*itv  cl *mnc»TioN  thu  »«oer*b«"  p»<»  t "««>«*) 


3 


20.  Abstract  (continued) 

and  dynamic  loading  and  the  free-field  (in-situ)  behavior  of  deep-based 
structures. 


% 


Results  show  that  the  primary  physical  parameter  in  structural  response 
is  the  lateral  confinement  of  the  test  specimen.  Tunnel  response  is  shown  to 
be  highly  sensitive  to  underconfinement,  although  relatively  insensitive  to 
moderate  overconfinement.  The  SRI  test  procedure  provides  slight  overconfine¬ 
ment  in  static  tests  designed  to  represent  the  uniaxial  strain  free-field 
conditions.  Therefore,  the  static  test  gives  a  good  estimate  of  tunnel 
closure  under  uniaxial  strain  conditions.  The  static  isotropic  (hydrostatic) 
loading  gives  smaller  (less  conservative)  tunnel  closures  than  uniaxial 
strain  loading. 


Tunnel  closures  in  dynamic  loading  tests  were  observed  to  be  much  smaller 
than  in  comparable  static  loading  tests.  Based  on  experimental  and  numerical 
results  the  reduction  in  dynamic  closures  are  most  likely  the  result  of  strain- 
rate  material  effect,  porewater  pressure  and  overconf inemept .  The  accuracy 
of  the  SRI  test  facility  in  indicating  the  response  of  actual  deep-based 
structures  in  a  dynamic  environment  cannot  be  determined  given  the  unresolved 
questions  of  material  behavior  and  confinement.  Because  of  this  uncertainty 
it  is  recommended  that  the  SRI  dynamic  tests  for  tunnel  closure  not  be  con¬ 
tinued.  However,  since  strain-rate  material  effects,  as  well  as  porewater 
pressure,  appear  to  be  significant,  it  is  recommended  that  static  and  dynamic 
tests  of  solid  specimens  of  site  material  and  rock  simulant  be  conducted  by 
SRI  to  provide  valuable  material  property  data. 
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SUMMARY 


S-l  BACKGROUND 

Test  facilities  have  been  developed  at  SRI  International  with 
the  objective  of  providing  inexpensive  small-scale  experimental 
techniques  capable. of  investigating  the  response  of  deep-based 
tunnel  structures  to  underground  loadings.  The  two  SRI  experi¬ 
mental  facilities  utilize  a  small  4  to  12  inch  diameter  cylindrical 
rock  specimen  into  which  is  cast,  perpendicular  to  the  axis  of  the 
cylinder,  a  5/8  to  2  inch  cylindrical  tunnel  liner,  as  shown  in 
Figure  S-l.  These  scale  models  of  tunnel/rock  structure  are  sub¬ 
jected  to  specified  (static  or  dynamic)  pressure  loads,  P^  and  P^ 
in  Figure  S-l,  on  the  ends  of  the  test  cylinder  and  on  its  lateral 
surface . 

This  report  presents  an  analysis  performed  in  support  of  the 
SRI  experimental  program  with  the  primary  objective  being  to 
evaluate  the  capability  of  the  test  apparatus  to  represent  the 
response  of  deep-based  tunnel  structures.  In  order  to  meet  this 
objective,  it  was  necessary  to  first  develop  an  understanding  of 
the  actual  experiments  and  to  then  isolate  the  most  significant 
aspects  of  the  SRI  test  procedure  as  a  representation  of  actual 
field  conditions. 

Of  the  various  modeling  questions  associated  with  the  SRI 
experiments  the  ones  which  were  believed  most  important,  and  which 
received  the  most  attention  in  the  present  investigation  were  the 
following : 

(1)  Given  the  applied  vertical  pressure,  what  is  the  effect 
of  lateral  confinement  and  how  is  this  effect  taken  into 
account  in  the  SRI  tests? 

(2)  What  is  the  sensitivity  of  the  material  variability 
on  the  test  results? 


(3)  What  is  the  significance  of  SRI  test  results  in  regard 
to  actual  deep-based  structural  response? 

Although  all  these  questions  pertain  to  the  SRI  tests,  the  first 
one  is  of  primary  consideration  in  evaluating  the  capability  of 
the  test  facility.  The  last  two  are,  in  addition,  important 
questions  which  must  be  dealt  with  in  order  to  design  deep-based 
structures  .* 

In  order  to  answer  these  questions  and  to  evaluate  the  SRI 
small-scale  test  facility,  a  comprehensive  analysis  program  was 
established  using  as  the  basis  the  CRT/NONSAP  finite  element  non¬ 
linear  analysis  code.  Several  finite  element  models  were  employed 
using  plane  strain  and  generalized  plane  strain  elastic-plastic 
representations  of  the  rock  simulant/tunnel  structure.  These 
computer  models  were  used  to  analyze  the  results  of  the  SRI  test 
facility  for  both  dynamic  and  static  loadings. 

S-2  DISCUSSION  OF  RESULTS 

Under  the  initial  program,  plane  strain  finite  element  anal¬ 
yses  of  certain  of  the  SRI  small-scale  tunnel  response  experiments 
were  performed.  The  general  objective  was  to  help  understand  and 
interpret  the  significance  of  some  of  the  observations  from  the 
nominally  uniaxial  strain  loading  tests.  In  particular,  substan¬ 
tially  smaller  tunnel  distortions  were  observed  in  dynamic - loaded 
tunnels  as  compared  to  static- loaded  tunnels. 

Figure  S-2  summarizes  the  most  important  results  of  the  pre¬ 
liminary  analyses.  Using  the  available  material  property  data 
(from  the  SRI  experiments) ,  the  plane  strain  analyses  essentially 
duplicated  the  tunnel  closure  versus  applied  stress  observed  by 
SRI  in  static  uniaxial  strain  test  SUX-103.  Then,  using  the  same 
properties  but  loading  the  tunnel  along  the  moderately  overconfined 
stress  path  which  occurred  in  SRI  dynamic  test  DUX- 77,  the  observed 
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Figure  S-2.  Preliminary  Analysis  Comparison  of  Experimental  and 
Theoretical  Crown  Closure  in  SRI  RMC  2C2. 
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closure  in  that  test  was  also  essentially  duplicated  by  numerical 
analyses.  In  both  cases,  the  results  may  have  been  fortuitous,  so 
additional  analyses  were  clearly  desirable  before  firm  conclusions 
could  be  drawn.  However,  these  results  are  self-consistent,  and 
they  do  suggest  that  the  differences  in  lateral  confinement  could 
be  responsible  for  the  differences  in  the  static-dynamic  results. 

As  a  corollary,  the  results  suggest  that  tunnel  response  may  be 
sensitive  to  the  degree  of  lateral  confinement.  This  could  be 
important  to  buried  structure  design  criteria. 

The  results  of  the  generalized  plane  strain  numerical  calcul¬ 
ations  demonstrate  the  usefulness  of  theoretical  analyses  for  inter¬ 
pretation  of  deep-based  tunnel  structures.  The  analytical  results 
are  in  good  agreement  with  experimental,  given  the  uncertainty  in 
material  behavior.  For  example,  Figures  S-3  and  S-4  compare  the 
numerical  and  experimental  crown  tunnel  closure  and  lateral  strain 
variation  with  lateral  confinement.  The  experimental  and  calculated 
results  show  similar  trends,  with  the  calculated  closure  being  con¬ 
sistently  less  than  the  experimental. 

These  differences  between  theoretical  and  experimental  results 
are  well  within  the  range  of  material  variability  typical  of  SRI 
RMG  2C2.  From  the  analyses  of  the  effect  of  material  behavior  on 
tunnel  response,  it  appears  that  (for  a  given  confinement)  material 
strain  rate  is  the  dominant  effect  in  reducing  tunnel  closure  in 
dynamic  tests  and  porewater  pressure  is  the  dominant  effect  in 
increasing  tunnel  closure  in  static  tests. 

Perhaps  the  most  significant  question  in  regard  to  the  SRI 
test  apparatus  is  its  capabilities  to  simulate  the  response  of  actual 
deep-based  tunnels.  To  evaluate  and  interpret  the  effectiveness 
of  the  SRI  experiment,  numerical  calculations  were  performed  to 
isolate  and  identify  discrepancies  that  could  occur  between  ideal 
free-field  tests  and  SRI  tests  performed  on  finite  samples.  Thus, 
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Figure  S-3.  Comparison  of  Experimental  and  Theoretical  Crown 
Closure  in  SRI  RMG  2C2  for  Uniaxial  Strain,  Over¬ 
confined  and  Underconfined  Loadings.  (Aluminum 
liner,  a/h  =  11.5.  is  lateral  pressure  required 

to  maintain  uniaxial  strain  as  defined  by  SRI.) 
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LATERAL  STRAIN  IPERCENT) 


several  analyses  were  performed  using  finite  element  grids  where 
the  boundaries  were  far  removed  from  the  tunnel.  These  analyses 
indicate  a  slightly  overconfined  condition  along  the  surface  of  the 
SRI  test  specimen,  as  compared  to  the  associated  free-field  surface. 

In  order  to  show  the  sensitivity  of  tunnel  closure  to  confine¬ 
ment,  the  results  of  the  present  study  are  re-plotted  in  Figure 
S-5  and  compared  with  four  definitions  of  neutral  confinement. 
Examining  Figure  S-S  shows  tunnel  closures  to  be  extremely  sensitive 
to  underconfinement  using  all  four  definitions  of  neutral  confine¬ 
ment.  These  results  also  show  that  the  definition  of  neutral  con¬ 
finement  can  significantly  influence  the  behavior  of  tunnel  closure 
(and  lateral  boundary)  response.  On  the  other  hand,  as  the  test 
specimen  is  overconfined,  the  tunnel  crown  closure  approaches  that 
for  the  ideal  uniaxial  strain  (roller)  boundary  state.  The  iso¬ 
tropic  (hydrostatic)  loading  condition  gives  smaller  (less  conser¬ 
vative)  tunnel  closures  than  uniaxial  strain  loading. 

These  analyses  support  the  general  validity  of  the  SRI  experi¬ 
ments  as  a  means  for  predicting  tunnel  response  of  an  idealized 
uniaxial  strain  test  under  a  static  loading  environment.  (However, 
some  of  the  results,  specifically  lateral  boundary  behavior  appear 
inconsistent.  For  example,  numerical  analyses  of  the  SRI  simulated 
uniaxial  strain  test  indicate  that  significant  lateral  bulging 
occurs  especially  at  high  pressure  levels.  This  leads  to  the 
obvious  question  as  to  how  the  SRI  test  appears  to  deviate  from  the 
ideal  uniaxial  strain  configuration,  yet  give  tunnel  closures  which 
are  essentially  equal.)  This  follows  from  the  fact  that  lateral 
boundary  behavior  has  only  a  secondary  effect  on  tunnel  closures 
when  overconfined .  Thus,  any  simulated  uniaxial  strain  test  that 
is  slightly  overconfined  will  essentially  duplicate  tunnel  response 
of  ideal  uniaxial  strain  tests. 
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S-3  RECOMMENDATIONS 


The  present  analyses  of  the  SRI  tests  support  the  capability 
of  the  uniaxial  strain  test  in  predicting  ideal  free-field  tunnel 
closure  for  static  loads.  These  tunnel  closure  results  are,  how¬ 
ever,  very  sensitive  to  lateral  (under)  confinement,  which  SRI 
controls  through  instrumentation  at  two  boundary  stations.  It  is 
recommended  that  this  lateral  confinement  be  more  carefully  con¬ 
trolled  by  instrumenting  the  specimens  at  several  stations  and 
controlling  maximum  lateral  deformation  to  insure  that  overcon¬ 
finement  will  always  be  present.  It  should  be  recognized  that  as 
overconfinement  increases,  vertical  tunnel  closure  decreases. 
Therefore,  in  order  for-  the  test  to  be  relatively  conservative,  it 
is  necessary  for  the  tests  to  be  only  slightly  overconfined. 

Since  lateral  confinement  in  the  SRI  dynamic  loading  test 
facility  is  difficult  if  not  impossible  to  control  and  since 
dynamical  material  property  data  is  not  avaiable  for  either  tuff 
or  rock  simulant,  it  is  not  possible  to  evaluate  the  capability 
of  the  SRI  test  to  simulate  far  field  dynamical  response.  Given 
these  uncertainties,  it  is  not  possible  to  support  additional 
dynamical  tests  for  tunnel  response  using  the  SRI  test  facility. 

However,  since  a  knowledge  of  rate-dependent  material  prop¬ 
erties  of  both  tuff  and  rock  simulant  is  of  fundamental  importance 
to  the  determination  of  deep-based  structural  response  to  dynamical 
loads,  it  is  recommended  that  the  SRI  facility  be  used  to  determine 
the  response  of  solid  specimens  to  dynamical  loads.  In  this  way 
a  data  base  can  be  established  which  will  be  essential  for  deter¬ 
mining  the  dynamical  material  properties  necessary  for  analyzing 
and  designing  dynamically  loaded  deep-based  structures. 
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PREFACE 


This  report  describes  a  two  part  investigation  performed  for 
the  Defense  Nuclear  Agency  under  Contracts  DNA001- 77-C-0248  and 
DNA001-78-C-0137  during  the  period  20  May  1977  through  15  October 
1978.  The  technical  monitors  for  this  investigation  were 
LTC  John  C.  Galloway  and  Dr.  Kent  L.  Goering. 

The  authors  wish  to  acknowledge  the  support  provided  by 
Kenneth  N.  Kreyenhagen  who  supervised  this  investigation  and 
the  valuable  discussions  with  Herbert  E.  Lindberg  and  Paul  E. 
Senseny  of  SRI  International. 
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Conversion  Factors  for  U.S.  Customary 
to  Metric  (SI)  Units  of  Measurement 


To  Convert  From 

To 

Multiply  By 

angstrom 

meters  (m) 

1.000  000  X  E  -10 

atmosphere  (normal) 

kilo  pascal  (kPa) 

1.013  25  X  E  +2 

bar 

kilo  pascal  (kPa) 

1.000  000  I  E  +2 

bam 

meter  (a  ) 

1.000  000  X  E  -28 

British  thermal  unit 
(thermochemical) 

joule  (J) 

1.054  350  X  E  +3 

calorie  (thermochemical) 

joule  (J) 

4.184  000 

cal  (thermochemical) /cm 

mega  joule /a*  (MJ/m*) 

4.184  000  X  E  -2 

curie 

glga  becquerel  (GBq)* 

3.700  000  A  E  +1 

degree  (angle) 

radian  (rad) 

1.745  329  X  E  -2 

degree  Fahrenheit 

degree  kelvin  (K)  * 

XK  -  (t*  f  +  *.59.67)/l.l 

electron  volt 

joule  (J)  | 

1.602  19  X  E  -19 

erg 

Joule  (J)  1 

1.000  000  X  E  -7 

erg/second 

ware  (w) 

1.000  000  X  E  -7 

foot 

meter  (m) 

3.048  000  X  E  - 1 

foot-pound- force 

Joule  (J) 

1.355  818 

gallon  (U.S.  liquid) 

meter*  (m*) 

3.785  412  X  E  -3 

inch 

meter  (m) 

2.540  000  X  E  -2 

Jerk 

joule  (J) 

1.000  000  X  E  49 

joule /kilogram  (J/kg> 

(radiation  dose  absorbed) 

Crey  (Gy)** 

1.000  000 

kilotons 

terajoules 

4.183 

kip  (1000  lbf) 

newton  (N) 

4.448  22 2  X  E  +3 

Up/inch2  (Its  1 J 

kilo  pascal  (kPa) 

6.894  757  X  E  +3 

ktap 

newtoo-eecond/m 

<N-s/aJ) 

1.000  000  X  E  +2 

micron 

meter  (a) 

1.000  000  X  E  -6 

mil 

meter  (a) 

2.540  000  X  E  -5 

mile  (International) 

meter  (m) 

1.609  344  X  E  +3 

ounce 

kilogram  (kg) 

2.834  952  X  E  -2 

pound- force  (lbf  avoirdupois) 

newton  (N) 

4.448  222 

pound- force  inch 

newton-meter  (N-tt) 

1.129  848  X  E  -1 

pound- force/ inch 

newton/meter  (N/m) 

1.751  268  X  E  +2 

pound-force/foot 

kilo  pascal  (kPa) 

4.788  026  X  E  -2 

pound- force/ Inch* 

kilo  pascal  (kPa) 

6.894  757 

pound-mass  (lbm  avoirdupois) 

kilogram  (kg) 

4.535  924  X  E  -1 

pound-mass- foot 

kilogram-meter* 

(moment  of  inertia) 

(kg/a3) 

1.601  646  X  E  +1 

rad  (radiation  dose  absorbed) 

Grey  (Gy)** 

1.000  000  X  E  -2 

roentgen 

coulomb /kilogram  (C/kg) 

2.579  760  X  E  -4 

shake 

second  (s) 

1.000  000  X  E  -8 

■  lug 

kilogram  (kg) 

1.459  390  IE  +1 

torr  (mi  Hg,  0*  C) 

kilo  pascal  (kPa) 

1.333  22  X  E  -1 

•The  becquerel  (Bq)  is  the  SI  unit  of  radioac t ivltv ;  1  Bq  •  l  event/s. 

••The  Cray  (Gy)  is  the  SI  unit  of  absorbed  radiation. 

A  more  complete  Listing  of  conversions  may  be  found  in  "Metric  Practice  Guide  E  380-74,'* 
American  Society  for  Testing  and  Materials. 
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SECTION  1 
INTRODUCTION 

1.1  REVIEW  OF  SRI  SMALL-SCALE  TEST  FACILITY. 

SRI  International  has  developed  [1-43*  laboratory- scale  test 
facilities  for  the  response  of  tunnel  liners  and  other  buried 
structures  (as  well  as  their  internal  components)  under  static  or 
dynamic  loading.  A  small-scale  test  apparatus  provides  a  cost- 
effective  way  of  screening  structural  concepts  and  conducting  para¬ 
metric  studies  of  buried  structural  response.  Figure  1.1  illus¬ 
trates  SRI's  small-scale  test  apparatus,  with  independent  appli¬ 
cation  of  vertical  and  lateral  pressure  under  both  static  and 
dynamic  load  environments.  The  initial  apparatus  was  designed 
[13  to  test  a  4-inch  diameter  cylinder  with  tunnel  models  of  approx¬ 
imately  5/8-inch  diameter,  as  shown  in  Figure  1.1.  Later,  a  larger 
apparatus  was  developed  C43  that  could  test  12-inch  diameter  cylin¬ 
ders  with  larger  diameter  tunnel  models. 

For  static  loading  (see  Figure  1.2),  the  pressures  on  the  top 
surface ,  PV>  and  on  the  sides  of  the  test  cylinder,  P^,  are  manually 
adjusted.  For  dynamic  loading,  explosives  are  used  indirectly  to 
generate  the  pressures  Pj.[  and  Py.  The  test  cylinder  "floats"  on 
a  lower  reservoir  of  oil  which  applies  (nominally)  equal  pressure 
Py  to  the  top  and  bottom  surfaces. 

The  apparatus  was  originally  designed  to  give  loading  rise 
times  of  the  order  of  50-100  psec.  Such  rapid  loading,  however, 
led  to  high  frequency  wave  reverberation.  Therefore,  the  apparatus 
was  modified  to  produce  relatively  gradual  dynamic  loads,  with 
rise  times  of  about  5  msec.  This  is  equivalent  to  a  rise  time 
of  perhaps  a  full  second  on  a  real  tunnel  liner.  Based  on  this 
modified  rise  time,  the  dynamic  tests  appear  to  be  quasi-static 
in  nature. 


*  Numbers  in  brackets  designate  references  at  the  end  of  the  report. 
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Isotropic 


Uniaxial 

Strain 


The  data  collected  from  a  dynamic  test  consist  of  pre-  and 
post-test  measurements*  of  the  crown- invert  and  springline  diameters, 
time  histories  of  the  vertical  and  lateral  pressures  and  time 
histories  of  hoop  strain  on  the  specimen's  lateral  surface.  In  a 
static  test,  the  vertical  and  lateral  pressures,  hoop  strain,  crown- 
invert  and  springline  diameters  were  recorded  during  each  increment 
(-500  psi)  of  loading. 

1.2  REVIEW  OF  SRI  TEST  PROCEDURE. 

The  primary  objective  of  the  SRI  test  fixture  is  to  simulate 
a  uniaxial  strain  condition  as  would  be  expected  in-situ.  It  was 
assumed  that  ideal  uniaxial  strain  of  the  finite  specimen  boundary 
was  a  reasonable  approximation.  In  static  tests,  the  test  fixture 
approximates  the  ideal  uniaxial  strain  lateral  boundary  using  the 
following  procedures: 

(1)  A  uniform  vertical  pressure  is  applied  hydraul¬ 
ically  to  the  cylindrical  test  specimen  while 
the  lateral  pressure  is  held  constant.  Vertical 
load  increments  are  applied  in  ~0.5  ksi  steps. 

(2)  For  each  vertical  load  increment  a  uniform 
lateral  pressure  is  then  applied  hydraulically 
until  the  lateral  boundary  indicates  zero 

lateral  motion.  Lateral  motion  is  approximated  by  the 
hoop  strain  on  the  lateral  boundary  as  measured  by 
the  average  of  two  strain  gauges  diametrically  opposed 
on  the  cylinder,  at  stations  one  tunnel  diameter  above 
the  mid-height. 


*  In  the  larger  test  apparatus, additional  stress  and  strain  gauges 
could  be  placed  at  various  locations  in  the  specimen.  Thus,  a 
dynamic  history  of  liner  deformations  could  be  made  based  on 
strain  gauges  placed  in  the  tunnel. 


Figure  1.3  illustrates  the  differences  that  can  exist  on  the 
lateral  boundary  for  an  ideal  uniaxial  strain  boundary  and  the 
SRI  static  test. 

In  the  dynamic  test,  the  vertical  and  lateral  loadings  are 
controlled  separately  but  with  a  fixed  ratio  using  low-density 
explosive  charges.  Unlike  the  static  test,  incremental  lateral 
loads  cannot  be  adjusted  to  maintain  a  zero  hoop  strain  measure¬ 
ment  on  the  lateral  boundary.  Instead,  the  uniaxial  strain 
(non-linear)  load  path  is  approximated  by  an  assumed  linear  load¬ 
ing  path  and  the  hoop  strains  are  recorded.  The  magnitudes  of 
these  lateral  strains  indicate  the  degree  of  overconfinement  or 
underconfinement  of  the  specific  test. 

1.3  REVIEW  OF  SRI  SMALL-SCALE  EXPERIMENTS. 

Using  the  small-scale  test  machine,  SRI  has  performed  [3,53 
tests  on  scale  models  of  MIGHTY  EPIC/DIABLO  HAWK  structures  in  a 
tuff  simulant,  SRI  RMG  2C2.  Both  dry  and  water- saturated  models 
were  tested  under  static  and  dynamic  conditions. 

These  test  specimens,  made  of  the  rock  or  tuff  matching 
grout,  provided  a  greater  specimen  homogeneity  than  found  in 
normal  geological  test  samples.  The  small-scale  test  cylinders 
were  4  inches  in  diameter  and  4  inches  high,  with  a  tunnel  bored 
through  the  cylinder  along  a  diameter  at  a  point  halfway  along 
its  height.  Each  of  the  several  tunnel  liner  models  tested  were 
inserted  into  the  tunnel  and  grouted  with  a  high  strength  gypsum 
cement . 

For  saturated  tests  in  which  porewater  drainage  was  permitted, 
eight  1/32-inch  holes  were  drilled  in  the  tunnel  liner  near  its 
ends  to  allow  the  seepage  of  porewater.  In  the  tests  where 
the  saturated  specimen  had  no  porewater  drain  holes,  a  pore¬ 
water  pressure  buildup  gave  anomalous  tunnel  deformations. 

Some  of  the  static  tests  had  drain  holes  only  at 
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one  end  of  the  liner  (creating  a  porewater  pressure  bulge  at  the 
undrained  end)  while  others  had  drain  holes  at  both  ends.  In 
dynamic  cases,  the  drainage  was  not  necessary  since  the  time  of 
porewater  migration  was  much  longer  than  the  time  frame  of  the 
test . 

Closure  data  of  SRI  RMG  2C2  is  shown  in  Figures  1.4  and  1.5. 

A  basic  trend  appears  in  the  data  showing  that  the  closures  for 
the  dynamic  tests  were  considerably  less  than  the  corresponding 
static  tests. 

1.4  REVIEW  OF  SRI  RMG  2C2  MATERIAL  DATA 

Since  the  rock  matching  grout,  SRI  RMG  2C2,  was  a  special 
material,  the  existing  constitutive  properties  were  not  readily 
available  and  had  to  be  determined  by  material  testing.  Initially, 
it  was  assumed  that  tunnel  response  would  not  influence  the  uni¬ 
axial  strain  behavior  of  the  SRI  test  specimens  and  that  these 
tests  could  be  used  to  determine  material  properties.  Using  the 
SRI  static  tests  of  tunnelled  specimens,  a  lateral  versus  vertical 
pressure  relation.  Figure  1.6,  was  generated  as  if  it  were  a  solid 
test  specimen  under  ideal  uniaxial  strain  conditions.  A  bilinear 
fit  based  on  an  elastic-perfectly  plastic  material  was  then  used 
to  approximate  the  SRI  data,  as  shown  in  Figure  1.6,  to  determine 
a  working  set  of  material  constants.  This  technique  is  recognized 
as  being  approximate  since  the  SRI  lateral  versus  vertical  pressure 
data  could  be  influenced  by  the  tunnel. 

Standard  material  tests  from  WES  and  Terra  Tek  were  also 
available.  These  tests  (triaxial  and  unconfined  compression  and 
uniaxial  strain)  were  conducted  on  solid  cylinders  of  SRI  RMG  2C2. 

Table  1.1  summarizes  the  interpretation  of  the  available  data, 
based  on  the  assumption  of  elastic-perfectly  plastic  behavior. 

The  differences  between  data  from  the  three  laboratories  are  not 
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Figure  1.5.  Tunnel  Closure  Versus  Applied  Pressure  for  Iso¬ 
tropic  Loading  of  SRI  RMG  2C2.  Liner:  6061-T0 
Aluminum,  a/h  =  11.5.  From  SRI  [5]. 
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Figure  1.6.  Lateral  Confining  Pressure  Versus  Vertical 
Pressure  Needed  to  Maintain  Uniaxial  Strain 
Conditions  in  SRI  RMG  2C2.  From  SRI  [5]. 


30 


Table  1.1.  Available  Material  Parameters  for  SRI  RMG  2C2. 


uniaxial  strain  test 
triaxial  failure  envelope 
test  specimen  contains  tunnel 


generally  large,  except  in  the  friction  angle,  which  is  much 
larger  in  the  WES  data.  This  observation,  combined  with  the  sub¬ 
stantially  greater  compaction  encountered  in  the  WES  data,  sug¬ 
gests  that  the  specimens  tested  by  WES  may  have  been  less  saturated 
than  the  samples  tested  at  SRI  and  Terra  Tek. 

In  summary,  these  tests  show  the  likelihood  of  a  moderate 
scatter  and  consequently  an  uncertainty  as  to  what  were  the  actual 
properties  in  the  SRI  RMG  2C2  when  tested  at  SRI. 

1.5  OBJECTIVES  OF  NUMERICAL  INVESTIGATION 

The  preliminary  SRI  results  indicated  some  sensitivity  to 
strain  rate  in  dynamic  tests  and  to  porewater  pressure  in  static 
tests.  A  major  reduction  in  tunnel  closure  (by  factors  of  10)  was 
observed  when  dynamic  tests  results  were  compared  to  similar  static 
tests.  Also,  porewater  pressure  buildup  gave  anomalous  tunnel 
deformations  in  static  tests  of  saturated  specimens. 

Clearly  these  observations  have  important  implications  with 
regard  to  both  the  SRI  apparatus  and  the  development  of  design 
data  for  buried  structures.  Also,  lateral  confinement  conditions 
are  important  but  rather  difficult  to  control.  (At  the  same  time 
variations  in  confinement  will  likely  occur  in  actual  geologic 
conditions  because  of  anisotropy,  geologic  variations,  proximity  of 
interfaces,  etc.) 

Changes  in  tunnel  closure  observed  for  dynamic  laboratory 
tests  versus  similar  static  tests  were  thought  to  be  due  to 
differences  in  lateral  confinement  and  to  material  property  vari¬ 
ations  from  strain-rate  and  porewater  pressure  effects. 

These  considerations  were  the  motivation  for  the  following 
objectives  of  the  present  investigation: 
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(1)  To  understand  and  interpret  the  results  of  the 
small-scale  SRI  tests. 

(2)  To  evaluate  the  capability  of  the  small-scale 
tests  to  simulate  deep-based  structural  response. 

(3)  To  determine  the  areas  of  greatest  sensitivity  in 
the  test,  namely  (a)  the  effect  of  lateral  confine¬ 
ment;  (b)  the  effect  of  material  property  variation 
and/or  definition. 

The  analyses  required  to  reach  these  objectives  are  detailed  in 
the  remainder  of  this  report.  The  CRT  version  of  the  NONSAP 
finite  element  computer  code  C6,7  3  was  used  to  perform  the  numer¬ 
ical  calculations. 

1.6  REPORT  ORGANIZATION 

The  next  two  chapters  give  the  numerical  models  used  to  per¬ 
form  the  present  analyses.  Chapter  2  describes  the  mathematical 
material  models  for  the  rock  simulant,  and  Chapter  3  gives  the 
finite  element  model  representations  of  the  uniaxial  strain  loading 
experiments.  Then,  Chapter  4  presents  the  preliminary  numerical 
analysis  of  the  SRI  experiments.  The  results  on  the  ideal  free- 
field  simulation,  the  effect  of  material  (including  rate)  behavior 
on  tunnel  response  and  the  sensitivity  of  tunnel  response  to 
lateral  confinement  are  given  in  Chapters  5,  6,  and  7,  respectively. 
Finally,  Chapter  8  presents  the  numerical  simulation  of  the  static 
uniaxial  strain  loading  experiments. 


SECTION  2 


MATHEMATICAL  MODELS  FOR  ROCK  SIMULANT,  SRI  RMG  2C2 

2.1  BACKGROUND 

In  the  initial  (preliminary)  phase  of  the  numerical  analyses, 
the  rock  simulant,  SRI  RMG  2C2,  was  idealized  with  an  elastic 
perfectly-plastic  material  model  using  the  Drucker -Prager  yield 
criterion  and  its  associated  flow  rule.  The  material  parameters  for 
the  model  were  determined  from  Py-P^  experimental  data  on  tunnelled 
cylinder  tests  provided  by  SRI.  The  material  constants  (cohesion 
and  friction  angle)  were  evaluated  assuming  that  the  tunnel  had  an 
insignificant  effect  on  the  boundary  response  (i.e.,  Py-P^  relation¬ 
ships  for  solid  cylindrical  specimens  correspond  to  Py-P^  relation¬ 
ships  of  tunnelled  specimens).  However,  the  preliminary  calcula¬ 
tions  (Section  4)  did  not  support  the  validity  of  this  assumption. 

That  is,  based  on  the  finite  sample  size  used  in  the  SRI  tests, 
lateral  boundary  behavior  is  influenced  by  the  presence  of  the 
tunnel.  Thus,  the  material  constituents  used  in  the  preliminary 
analysis  did  not  necessarily  reflect  the  characteristics  of  the 
test  material. 

For  subsequent  phases  of  the  analyses  (Sections  5  to  8) ,  material 
constituents  and  the  mathematical  description  thereof  were  evaluated 
from  test  data  corresponding  to  experiments  performed  on  solid 
cylinders.  Triaxial  compression  and  uniaxial  strain  test  data 
provided  by  WES  and  Terra  Tek  was  used  to  formulate  material  consti¬ 
tuent  models  and  to  evaluate  their  material  parameters. 

2.2  DRUCKER- PRAGER  MODEL 

The  Drucker-Prager  model  based  on  the  Coulomb  rule  for  slipping, 
has  a  conical  yield  surface  in  principal  stress  space.  The  yield 
function  resembles  von  Mises  in  the  ir-plane  but  with  a  pressure  term 
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added  to  account  for  the  Coulomb  effects  C8i 

F  =  aJi  +  /jT  -  k  (2.1) 

where  a  and  k  are  positive  constants  and  Ji  and  are  the  first 

stress  invariant  and  second  deviatoric  stress  invariant,  respectively. 

Drucker  and  Frager  proposed  C9l  this  model  for  limit  design,  so  they 

considered  rigid-plastic  strain  conditions,  i.e.,  the  only  non-zero 

strains  are  plastic  strains.  For  this  case,  they  showed  that  this 

yield  function  fits  the  Mohr-Coulomb  yield  condition  for  the  case 

P  F  P 

of  plastic  plane  strain,  i.e.,  e,^  =  e^3  =  e23  =  if  the  a  anc* 
k  values  were  taken  as : 


sin<?> 


and  k  =  _3_c_cos$_ 


(2.2) 


'9  +  3sin‘ 


'9  +  3sin' 


where  4>  and  c  are  properties  of  friction  angle  and  cohesion, 
respectively.  If  triaxial  test  results  are  used  to  find  4>  and  c, 
then  the  plastic  plane  strain  condition  is  no  longer  valid  and  a  new 
Drucker-Prager  condition  is'  derived.  This  will  be  designated  as  a 
triaxial  Drucker-Prager  condition  to  clarify  it  from  the  original 
Drucker-Prager  plane  strain  condition.  Using  triaxial  test  results, 
the  a  and  k  terms  take  the  form: 


2sin<{> 

/J- (3-sini}>) 


6  c  cost}) 

/3~  (3-sini))) 


(2.3) 


In  comparing  several  of  these  surfaces  at  their  Ti-plane  C  8  J  inter 
sections,  see  Figure  2.1,  it  should  be  noted  that  the  triaxial  Drucker 
Prager  circumscribes  the  Mohr-Coulomb  surface,  while  the  plane 
strain  Drucker-Prager  inscribes  it.  The  plane  strain  Drucker-Prager 
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surface  matches  the  Mohr-Coulomb  surface  at  the  rigid  plastic  plane 
strain  points,  while  the  triaxial  Drucker- Prager  model  matches  the 
Mohr-Coulomb  at  the  triaxial  compression  locations. 

Because  of  the  nature  of  the  SRI  tests  where  some  plastic  flow 
occurs,  and  triaxial  compression  governs,  the  triaxial  Drucker- 
Prager  model  was  used.  (The  NONSAP  Drucker- Prager  has  the  triaxial 
formulation 'built  into  it,  although,  the  plane  strain  formulation 
is  easy  to  incorporate.) 

The  validity  of  the  Drucker-Prager  model  was  verified  under 
uniaxial  strain  conditions  by  applying  stress  and  roller  boundary 
conditions.  The  NONSAP  results  were  then  compared  against  the 
analytical  solution  as  shown  in  Figure  2.2. 

This  Drucker-Prager  model  representation  of  the  SRI  RMG  2C2 
material  was  used  in  the  preliminary  numerical  simulation  of  the 
SRI  uniaxial  strain  experiments  given  in  Section  4. 

2.3  VARIABLE  FRICTION  ANGLE,  STRAIN  HARDENING  DRUCKER-PRAGER 

MATERIAL  MODEL 

When  typical  experimental  data  for  SRI  RMG  2C2  is  compared  with 
the  modeling  capabilities  of  either  a  per fectly-plastic  Drucker- 
Prager  or  Mohr-Coulomb  model,  the  following  deficiencies  become 
apparent : 

(1)  Figure  2.3  indicates  typical  experimental  values  of 

r 2  * 

octahedral  shear  stress,  /  yJ2 ,  as  a  function  of  mean  pressure 
for  both  initial  yield  and  failure.  This  data  was  obtained 
from  triaxial  compression  tests.  Experimental  data  shows 
that  both  the  initial  yield  and  failure  surfaces  vary  non- 
linearly  with  respect  to  mean  pressure.  Thus,  Drucker- 
Prager  or  Mohr-Coulomb  model  at  best  can  only  fit  the  data 
for  a  small  range  of  mean  pressure. 
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Figure  2.2.  Preliminary  Analysis  Comparison  of  Analytical  and 
Numerical  Results  for  the  Drucker-Prager  Model. 


Figure  2.3.  Comparison  of  Typical  Experimental  Failure  and  Yield 

Surface  Data  with  a  Perfectlv-Plastic  Drucker-Praeer  Model 


Figure  2.4.  Comparison  of  Typical  Unconfined  Stress-Strain  Data  with 
Perfectlv-Plastic  Drucker-Praeer  and  Mnlu — Foul  nmh  Models. 
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(2)  A  typical  unconfined  compression  test  for  SRI  RMG  2C2 
is  shown  in  Figure  2.4.  After  the  onset  of  plastic 
flow,  the  data  illustrates  the  strain -hardening  char¬ 
acteristics  of  SRI  RMG  2C2.  Strain-hardening  behavior 
of  SRI  RMG  2C2  is  also  indicated  in  the  expansion  from 
the  virgin  yield  surface  to  the  failure  surface  (See 
Figure  2.3).  The  disadvantage  of  a  perfectly-plastic 
model  is  that  it  does  not  allow  for  yield  surface  to 
change  either  in  size  or  shape. 

(3)  The  applied  stress -confinement  stress  ratios  for  a 
uniaxial  strain  test  of  SRI  RMG  2C2  is  shown  in  Figure  2.5. 
The  test  data  indicates  that  after  yielding,  the  applied- 
confinement  pressure  relation  varies  non- 1 inearly  .  Per¬ 
fectly-plastic  Drucker-Prager  and  Mohr-Coulomb  material 
models  do  not  predict  this  non-linear  behavior,  in  fact 
these  models  lead  to  bilinear  applied-confinement  stress 
relationships . 

(4)  The  numerical  simulation  in  the  inelastic  regime  of  a 
uniaxial  strain  test  by  the  application  of  confinement 
pressure  to  an  unconfined  specimen  is  not  possible  using 
perfectly-plastic  material  models.  When  yielding  occurs 
in  a  geometrically  unconfined  specimen  (uniaxial  simu¬ 
lation  provided  through  the  application  of  lateral  con¬ 
finement  pressure) ,  the  resulting  system  stiffness  is 
semi -definite ,  so  that  complete  collapse  results.  Figure 
2.6  shows  the  applied-confinement  pressure  that  exists  in  a 
state  of  uniaxial  strain  assuming  a  perfectly-plastic 
material  model.  Any  value  of  confinement  stress  that  lies 
slightly  to  the  right  of  the  uniaxial  strain  results  will 
produce  collapse  of  a  simulated  uniaxial  strain  test.  On 
the  other  hand,  applied-confinement  stress  relations  that 
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Figure  2.6.  Uniaxial  Strain,  Applied  Stress-Confinement  Stress  for  Perfectly 
Plastic  Drucker-Prager  and  Mohr-Coulomb  Material  Models 


lie  to  the  left  cause  the  geometrically  unconfined  specimen 
to  behave  elastically.  This  will  cause  the  simulated  uniax¬ 
ial  strain  to  deviate  substantially  from  the  uniaxial  strain 
state.  The  nature  of  the  SRI  test  implies  that  the  SRI  RMG 
2C2  must  strain-harden.  Otherwise,  the  SRI  static  load 
controlled  test  for  a  per  feet ly-plastic  or  strain-softening 
material  would  result  in  collapse. 

The  above-mentioned  inadequacies  of  a  perfect -plastic  Drucker- 
Prager  or  Mohr-Coulomb  model  for  the  idealization  of  SRI  RMG  2C2 
led  to  the  development  of  an  improved  plastic  material  model  with  the 
following  capabilities: 

(1)  Both  the  yield  surface  and  failure  surface  indicated  by 
the  experimental  data  can  vary  with  increasing  pressure 
in  a  non-linear  manner. 

(2)  The  material  model  has  work  dependent  variable  isotropic 
hardening  which  gives  improved  simulation  of  experimental 
data  and  permits  the  numerical  modeling  of  a  geometrically 
unconfined  simulated  uniaxial  test. 

The  details  involved  with  NONSAP  implementation  of  the  variable 
Drucker- Prager  model  and  a  comparison  of  its  fit  with  experimental 
data  is  given  in  Appendix  A.  This  model  is  used  in  Sections  5  to  8. 

2.4  MODIFIED  MOHR-COULOMB  MATERIAL  MODEL 

SRI  has  also  performed  theoretical  NONSAP  analyses  to  help  in 
their  evaluation  and  interpretation  of  the  experimental  results. 

Their  representation  of  the  material  constitutive  behavior  of  the 
rock  medium,  SRI  RMG  2C2,  consisted  of  an  elast ic -perfectly  plastic 
material  model  using  the  Mohr-Coulomb  yield  criterion  with  an  assoc¬ 
iated  plastic  flow  rule.  The  Mohr-Coulomb  material  subroutine  i  10  1 
was  added  to  the  NONSAP  code  by  SRI  for  the  specific  purpose  of 
performing  these  calculations.  SRI  completed  theoretical  calcula- 
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tions  for  isotropic ?  loading,  which  compared  reasonably  well  with 
static  experimental  results.  However,  SRI  had  difficulties  when 
they  attempted  to  apply  NONSAP,  specifically  their  Mohr-Coulomb 
material  model,  to  numerically  simulate  their  uniaxial  strain  tests. 

The  Mohr-Coulomb  material  model  was  modified  in  the  current 
program  to  alleviate  the  difficulties  encountered  by  SRI.  In  this 
modification. the  corners  of  the  traditional  Mohr-Coulomb  yield 
criterion  were  rounded,  so  that  the  yield  criterion  has  continuous 
normal  derivatives,  and  hence  unique  plastic-flow  conditions  for  all 
states  of  stress  (see  Ref.  r8l).  The  problem  encountered  by  SRI 
was  a  consequence  of  non-unique  plastic  flow  that  occurred  in  the 
corners  of  the  conventional  Mohr-Coulomb  yield  surface.  To  see  this 
effect,  consider  the  state  of  stress  under  uniaxial  strain  conditions 
(oj  =  Py,  Op  =  o3  =  P^) .  When  the  initial  yielding  occurs,  the  stress 
state  lies  at  a  corner  of  the  traditional  Mohr-Coulomb  surface.  Thus, 
any  round-off  error  introduced  by  a  computer  will  produce  any  arbitrary 
value  for  the  gradient,  introducing  asymmetry  and  numerical  instability 
in  the  problem.  The  modified  Mohr-Coulomb  criterion  eliminates  this 
problem  by  the  rounding  of  the  corners.  The  details  of  this  model 
are  described  in  Appendix  B. 
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SECTION  3 


FINITE  ELEMENT  MODELS  OF  UNIAXIAL 
STRAIN  LOADING  EXPERIMENTS 

The  intent  of  the  SRI  uniaxial  strain  loading  test  is  to 
represent  the  response  of  buried  structures  using  small-scale 
finite  specimens.  The  validity  of  the  SRI  test  is  in  a  sense 
dependent  on  the  assumption  that  the  finite  boundary  (approximately 
6  tunnel  diameters)  behavior  in  the  free-field  is  essentially 
equivalent  to  the  boundary  behavior  of  the  test  specimen  (see 
Figure  3.1).  The  SRI  test  was  conducted  under  the  assumption  that 
for  practical  purposes  a  uniaxial  strain  condition  at  the  finite 
boundary  is  equivalent  to  the  free-field  condition  (which  implies 
that  the  tunnel  response  does  not  influence  the  behavior  on  the 
finite  boundary).  Hence,  uniform  lateral  stress  is  applied  and  a 
uniaxial  strain  condition  is  prescribed  (actually  hoop  strain)  at 
a  point  along  the  boundary  (see  Figure  1.3). 

In  the  initial  (preliminary)  phase  of  the  analysis,  numerical 
studies  were  performed  (see  Section  4)  with  a  plane  strain  finite 
element  representation  of  the  SRI  test.  The  steps  employed  in  this 
analysis  are  described  as  follows.  First,  the  results  of  the  SRI 
simulated  uniaxial  strain  tests  (applied  pressure-confinement  pres¬ 
sures,  PV'PH)  were  used  to  evaluate  the  material  parameters  for  a 
perfectly-plastic  Drucker-Prager  model  as  discussed  in  Section  2.2. 
(The  validity  of  this  assumption  requires  the  tunnel  to  have  an 
insignificant  effect  on  the  response  at  the  boundary  of  SRI  test 
specimen.  That  is,  the  Py-P^  relation  for  a  uniaxial  strain  test 
of  a  solid  specimen  would  be  essentially  equivalent  to  that  with  a 
tunnel.)  Then,  having  established  the  material  properties,  the 
plane  strain  numerical  simulation  of  the  SRI  test  was  performed  by 
applying  the  experimental  py'pH  pressure  relation  as  stress  boundary 
conditions.  (With  this  stress  condition  no  lateral  strain  reference 
location  was  used.) 
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Figure  3.2  illustrates  the  plane  strain  modeling  capabilities 
and  limitations  with  regard  to  the  SRI  test.  The  plane  strain 
simulation  can  be  visualized  as  a  rectangular  slab  with  a  lateral 
hole  through  the  center,  with  uniform  stresses  applied  to  the  top 
and  lateral  surfaces.  Because  of  plane  strain  assumption,  the 
resultant  force  in  the  plane  strain  direction  may  not  be  equal  to 
the  resultant  lateral  force.  This  is  not  necessarily  consisent 
with  the  SRI  test  where  both  are  essentially  equal.  This  incon¬ 
sistency  can  be  resolved  by  removing  the  plane  strain  constraint 
(i.e.,  f  0) ,  as  shown  in  Figure  3.2 

Also,  the  sensitivity  of  the  preliminary  analysis  to  lateral 
confinement  led  to  the  development  of  an  alternate  procedure  for 
applying  the  lateral  constraint.  The  procedure  used  in  the 
numerical  approach  was  identical*  to  the  experimental  procedure. 

That  is,  lateral  pressures  were  adjusted  until  zero  lateral  dis¬ 
placement  was  obtained  at  a  point  on  the  lateral  boundary  cor¬ 
responding  to  the  location  of  the  SRI  reference  strain  gauges. 

Thus,  in  the  subsequent  phase  of  the  analysis,  numerical 
simulations  of  the  SRI  tests  were  performed  using  a  generalized 
plane  strain  finite  element  model  with  material  parameters  evaluated 
independently  of  the  SRI  test  data.  The  analysis  incorporates  two 
basic  improvements  over  the  preliminary  analysis:  (a)  the  zero 
lateral  strain  location  on  the  lateral  boundary  matches  the  experi¬ 
mental  location  and  (b)  the  plane  strain  constraint  is  relaxed, 
i.e.,  f  0,  which  allows  the  out-of-plane  force  to  be  specified 
(in  this  case,  equal  to  the  resultant  lateral  force). 

*The  two  hoop  strain  gauges  located  on  the  lateral  boundary  of  the 
SRI  test  specimen  were  intended  to  provide  a  measure  of  the  lateral 
displacement.  However,  due  to  asymmetry,  one  gauge  gave  a  positive 
strain  while  the  other  gave  a  negative  strain  for  an  average  strain 
of  zero.  This  constraint  does  not  necessarily  imply  that  the 
lateral  displacement  at  this  location  is  zero. 


SECTION  4 


PRELIMINARY  NUMERICAL  ANALYSIS 
OF  SRI  EXPERIMENTS 

Under  the  first  phase  of  the  program,  preliminary  analyses 
were  made  in  order  to  develop  a  numerical  model  of  the  SRI  test 
and  to  find  the  sensitive  factors  governing  test  data,  Figure  4.1. 

The  preliminary  model  used  a  plane  strain  finite  element 
representation  of  the  SRI  test  specimen  (see  Section  3) .  At  the 
lateral  boundaries  either  a  roller  boundary  condition  could  be 
applied  (an  ideal  uniaxial  strain  situation  in  a  solid  specimen) 
or  a  uniform  pressure  could  be  applied  based  on  the  SRI  test  (see 
Figure  4.2).  The  material  model  used  was  the  Drucker-Prager  model 
described  in  Section  2.2  and  shown  in  Figure  4.2. 

Analyses  were  performed  with  the  ideal  uniaxial  strain  (roller) 
conditions  on  the  lateral  boundary  and  compared  with  a  uniform 
stress  boundary  condition,  as  shown  in  Figure  4.3.  Although  a 
roller  boundary  condition  gives  a  state  of  uniaxial  strain  on  the 
sample,  the  presence  of  a  tunnel  alters  the  lateral  pressure  con¬ 
dition.  Figure  4.4  shows  the  results  of  the  analyses  performed 
to  evaluate  the  effects  that  tunnel  size  has  on  closure. 

Several  sets  of  material  properties  representing  both  dry  and 
water  saturated  SRI  RMG  2C2  were  simulated.  The  material  data  were 
obtained  from  P^-P^  curves  from  SRI's  static  tests.  Because  the 
SRI  dynamic  tests  were  essentially  quasi- static ,  the  tests  were 
initially  analyzed  by  duplicating  the  load  paths  in  a  static  solu¬ 
tion.  The  load  paths  were  taken  from  the  SRI  dynamic  test  data 
and  applied  to  the  plane  strain  grid,  so  that  the  amount  of  stress 
confinement  occurring  in  the  SRI  dynamic  test  is  duplicated.  In 
addition,  a  straight  line  approximation  to  the  load  path  was  used 
as  a  check  on  path  dependence.  These  results  are  shown  in  Figure 
4.5.  Since  the  dynamic  test  was  more  overconfined  than  the  static 
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Figure  4.2.  CRT/NONSAP  Material  Model  Fit  to  SRI  Test  Data  on 
Dry  Sample  (SUX-103)  . 
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Figure  4.3.  Numerical  Simulation  of  Static  Test  SUX-103,  using 
Roller  and  Pressure  Boundary  Conditions. 
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test,  as  shown  in  Figure  4.6,  the  results  should  be  essentially 
path  independent,  as  indicated  in  Section  7. 

The  preliminary  analysis  exposed  several  trends  in  the  SRI 
test  apparatus  (see  Figure  4.7).  These  initial  conclusions  were 
then  examined  more  carefully  in  the  second  phase  of  the  program. 

The  first  conclusion  reached  was  that  measurement  of  non¬ 
linear  material  properties  and  of  pressure  ratios  required  to 
maintain  uniaxial  strain  by  using  test  cylinders  containing  tun¬ 
nels  gives  misleading  information.  Thus,  tunnelled  specimens 
should  not  be  used  to  simultaneously  define  material  properties 
and  analyze  the  SRI  test.  Clearly,  solid  test  specimens  should 
have  been  used  in  the  SRI  test  apparatus  to  first  define  the 
material  behavior.  As  a  consequence  it  appears  that  the  ratio  of 
test  cylinder  diameter  to  tunnel  diameter  effects  the  closure, 
especially  after  the  onset  of  material  plasticity. 

Static  tests  should  be  emphasized,  at  least  for  now.  They 
are  much  less  expensive,  yield  more  complete  data,  are  more 
reliable,  and  are  inherently  easier  to  specify  and  control.  As 
such,  static  tests  are  a  good  way  to  ferret  out  such  things  as 
the  possible  sensitivity  to  lateral  confinement. 

Dynamic  loading  has  been  harder  to  control,  and  most  of  the 
dynamic  tests  have  been  overconfined  (as  compared  to  the  nominally 
uniaxial  strain  loading  in  the  static  tests)  . 

Figure  4.7  summarizes  the  most  important  results  of  the  pre¬ 
liminary  analyses.  Using  the  available  material  property  data, 
the  plane  strain  finite  element  analyses  essentially  duplicated 
the  tunnel  closure  versus  applied  stress  observed  by  SRI  in  static 
uniaxial  strain  test  SUX-103.  Then,  using  the  same  properties 
but  loading  the  tunnel  along  the  moderately  overconfined  stress 
path  which  occurred  in  SRI  dynamic  test  DUX-77,  the  observed 
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Figure  4.7.  Preliminary  Analysis  Comparison  of  experimental  and 
Theoretical  Crown  Closure  in  SRI  RMG  2C2. 
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closure  in  that  test  was  also  essentially  duplicated  by  numerical 
analyses.  In  both  cases,  the  results  may  have  been  fortuitous,  so 
additional  analyses  are  clearly  desirable  before  firm  conclusions 
are  drawn.  However,  these  results  are  self-consistent,  and  they 
do  suggest  that  the  differences  in  lateral  confinement  could  be 
responsible  for  the  differences  in  the  static-dynamic  results. 

As  a  corollary,  the  results  suggest  that  tunnel  response  may  be 
sensitive  to  the  degree  of  lateral  confinement.  This  could  be 
important  to  buried  structure  design  criteria. 


SECTION  5 


IDEAL  FREE  -  FIELD  SIMULATION 

Perhaps  the  most  significant  question  in  regard  to  the  SRI 
test  apparatus  is  its  capabilities  to  simulate  the  response  of  actual 
deep-based  tunnels.  To  evaluate  and  interpret  the  effectiveness  of 
the  SRI  experiment,  numerical  calculations  were  performed  to  isolate 
and  identify  discrepancies  that  could  occur  between  ideal  free-field 
tests  and  SRI  tests  performed  on  finite  samples.  Thus,  several 
analyses  were  performed  using  grids  where  the  boundaries  were  far 
removed  from  the  tunnel  (ratio  of  tunnel  diameter  to  cylinder  di¬ 
ameter  of  18)  and  the  improved  material  model  given  in  Appendix  A. 

In  the  SRI  uniaxial  strain  test,  uniform  stress  is  applied  to  the 
lateral  boundary  at  ~6a  so  that  nearly  zero  strain  occurs  (2a  above 
the  mid-height).  Using  the  results  of  the  far-field  boundary  anal¬ 
ysis,  the  stress  and  strain  distribution  along  the  6a  lateral 
boundary  can  be  compared  with  the  SRI  test. 

A  comparison  of  the  lateral  boundary  conditions  for  the  SRI 
test  with  free-field  conditions  is  a  key  consideration,  since  these 
boundary  conditions  have  a  very  strong  influence  on  the  specimen 
response.  Nevertheless,  the  fundamental  question  is  how  do  the 
tunnel  closures  for  the  SRI  test  compare  with  free-field?  Is  the 
SRI  test  overconfined  (non-conservative)  or  underconfined  (con¬ 
servative)  ? 

The  numerical  results  obtained  for  tunnel  crown  and  springlinc 
closures  for  the  simulated  free-field  condition  are  shown  in  Figures 
5.1  and  5.2,  respectively.  Three  different  assumptions  regarding 
the  nature  of  the  lateral  boundary  conditions  in  the  far-field 
(ratio  of  tunnel  diameter  to  cylinder  diameter  of  18)  were  used: 

(1)  Plane  strain  analysis  assuming  uniaxial  strain  (roller) 
condition  on  the  far-field  lateral  boundary. 
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(P*  is  the  confinement  pressure  which  maintains  uniaxial  strain 
cohdition  for  a  tunnelless  analysis.) 


Figure  3.2.  Comparison  of  Tunnel  Springline  Closures  for  Plain 
Strain  and  Generalized  Plane  Strain  Analyses  with 
Mound. lry  Cond i t  ions  Appl  iod  at  18a. 


(2)  Generalized  plane  strain  analysis  with  a  uniaxial  strain 
(roller)  condition  on  the  far-field  lateral  boundary 
where  the  average  out-of-plane  stress  is  equal  in  mag¬ 
nitude  to  the  ideal  free-field  lateral  stress. 

(3)  Generalized  plane  strain  analysis  with  uniform  lateral 
stress  applied  to  the  far-field  lateral  boundary  where 
both  the  lateral  stress  and  the  average  out-of-plane 
stress  are  equal  to  the  ideal  free-field  lateral  stress. 

The  closure  results  are  nearly  identical  for  the  two  cases  in 
which  the  lateral  boundary  is  confined  with  a  roller  boundary;  how¬ 
ever,  closure  increases  when  the  free-field  lateral  stress  boundary 
is  used.  This  clearly  suggests  that  a  relatively  small  change  to 
underconfinement  causes  rather  large  changes  in  tunnel  closures. 

The  non-uniform  stress  and  strain  distributions  at  the  6a 
boundary  obtained  for  the  uniaxial  strain  boundary  condition  at  18a 
are  shown  in  Figures  5.3  and  5.4,  respectively.  Figure  5.4  indi¬ 
cates  that  in  the  SRI  tests  an  improved  idealization  of  the  free- 
field  test  could  result  by  zeroing  the  strain  along  the  lateral 
boundary  at  a  location  5a  above  the  mid-height,  instead  of  2a. 

Thus,  the  SRI  test  procedure  is  essentially  overconfined  and  non¬ 
conservative  . 
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SECTION  6 


EFFECT  OF  MATERIAL  BEHAVIOR 
ON  TUNNEL  RESPONSE 

This  section  presents  the  numerical  analyses  of  tunnel  response 
to  variations  in  material  behavior.  Variability  in  material  char¬ 
acteristics  for  RMG  2C2  (material  used  in  the  SRI  static  and 
dynamic  tests)  appears  to  be  important,  based  on  deviations  in 
tunnel  response  that  occurred  in  "identical"  static  and  dynamic 
tests.  The  numerical  analyses  show  the  sensitivity  of  tunnel 
response  to  material  strength  (i.e.,  SRI  test  sensitivity  to  RMG 
2C2  strength).  In  addition,  they  provide  data  on  the  variations 
in  tunnel  response  which  are  likely  to  occur  from  differences  in 
RMG  2C2  and  actual  tuff  material  strengths. 

6.1  EFFECT  OF  MATERIAL  PROPERTY  VARIATIONS 

Triaxial  compression  tests  on  material  samples  taken  from  the 
MIGHTY  EPIC/DIABLO  HAWK  site  indicate  considerable  variability. 
Figure  6.1  illustrates  the  scatter  in  the  failure  surface  data 
corresponding  to  51  triaxial  compression  tests  of  MIGHTY  EPIC  tuff. 
The  material  selected  for  the  majority  of  the  testing  performed  at 
SRI  was  RMG  2C2,  a  tuff  simulant.  RMG  2C2  has  a  material  behavior 
which  reasonably  approximates  MIGHTY  EPIC  tuff  for  a  wide  range  of 
pressure.  Nevertheless,  the  geological  heterogeneity  in  MIGHTY 
EPIC  tuff  will  inevitably  be  greater  than  RMG  2C2,  hence  the  varia¬ 
tion  in  tunnel  closures  in  the  SRI  tests  are  likely  to  be  less  than 
in  actual  site  tests. 

Figure  6.2  shows  the  scatter  in  failure  surface  data  corres¬ 
ponding  to  five  separate  tests  performed  on  RMG  2C2  by  WES.  Also 
indicated  in  this  figure  are  the  three  sets  of  material  properties 
used  for  the  present  sensitivity  study.  The  variability  in  the 
assumed  material  properties  exceeds  the  scatter  in  the  WES  experi¬ 
mental  data.  However,  the  WES  tests  were  performed  on  specimens 


6  1 


■  IIH 


TERIW  TEK  DflTH 

FRItURE  ENVELOPES  FOR  TUFF  SIHULANT3 


3 it iv itt  effects  caused  ar  material  variations 

/SHIP 


with  unknown  water  content  under  static  conditions,  whereas  in  the 
SRI  tests,  the  water  content  varied  from  dry  to  saturated  under 
static  and  dynamic  conditions.  The  three  sets  of  material  prop¬ 
erties  reflect  assumed  nominal  (dry) ,  strong  (strain  rate  effect) 
and  weak  (porewater  effect)  behavior. 

The  numerical  analyses  regarding  the  effects  of  material 
variability  on  tunnel  closure  were  performed  using  a  plane  strain 
model  with  rigid  confinement  on  the  lateral  boundary  at  r  =  6.4a, 
the  numerical  simulation  of  an  ideal  uniaxial  strain  boundary  at 
6.4a.  The  results  of  this  study  are  summarized  in  Figures  6.3 
through  6.5.  In  Figure  6.3,  the  appl ied- confinement  pressures  for 
the  material  parameter  study  are  compared  with  the  SRI  experimental 
data  (the  analytical  lateral  pressure  represents  the  average 
lateral  stress).  Figures  6.4  and  6.5  compare  tunnel  crown  and  spring¬ 
line  closures.  The  numerical  predictions  agree  quite  well  with 
the  experimental  measurements  for  both  confinement  and  closure. 

The  calculated  results  not  only  match  the  experimental,  but 
also  the  band  in  experimental  data  matches  reasonably  well  with 
the  one  for  the  assumed  material  variations.  Actually,  the  scat¬ 
ter  in  the  experimental  data  is  more  than  likely  a  consequence  of 
a  small  variation  in  lateral  confinement,  rather  than  a  material 
variation.  The  variability  in  material  strength  used  in  this 
study  represents  a  reasonable  set  of  bounds  to  tuff  found  at  the 
MIGHTY  HPIC  site.  In  summary,  these  analyses  indicate  that  the 
scatter  in  SRI  test  data  caused  by  material  variations  and/or 
slight  variation  in  lateral  confinement  are  likely  to  be  /ra*? 
than  those  observed  in-situ. 

6.2  RATH  EFFHCTS 

SRI  test  data  for  tunnel  closure  indicates  that  major  dif¬ 
ferences  occur  under  quasi-static  and  dynamic  loading  environments. 
Figure  4.1  shows  typical  static  and  dynamic  test  data  for  uniaxial 
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strain  loading  of  SRI  RMG  2C2.  This  figure  clearly  indi¬ 
cates  the  relative  reduction  in  tunnel  distortion  for  dynamically 
loaded  specimens.  The  loading  in  the  SRI  simulated  dynamic  uniaxial 
strain  test  is  difficult  to  control,  and  most  dynamic  uniaxial 
strain  tests  were  overconfined.  Overconfinement  results  in  less 
closure  and  may  be  responsible  for  at  least  a  portion  of  the  dif¬ 
ferences  in  closures  between  static  and  dynamic  uniaxial  strain 
tests.  Nevertheless,  substantial  differences  in  tunnel  distortions 
occurred  for  static  and  dynamic  testing  under  isotropic  loading  (see 
Figure  1.5),  where  dynamic  loading  could  be  controlled.  This 
implies  that  the  differences  in  static  and  dynamic  uniaxial  strain 
tests  are  not  merely  the  product  of  overconfined  dynamic  tests. 


Typical  rise  times  in  the  SRI  dynamic  tests  are  on  the  order 
of  2  msec.  Based  on  the  material  properties  and  the  dimensions 
of  the  test  specimens  (4- inch  diameter)  the  rise  times  for  the 
dynamic  test  loading  are  very  slow  compared  to  wave  speeds.  Hence, 
inertia  effects  should  not  be  significant.  However,  the  blast 
loading  in  the  dynamic  tests  produces  strain  rates  several  orders 
of  magnitudes  greater  than  static  tests.  These  high  strain  rates 
may  effectively  strengthen  the  material  and  thereby  cause  less 
closure  in  dynamic  tests. 


The  procedure  employed  by  SRI  for  the  measurement  of  tunnel 
distortion  varies  under  static  and  dynamic  loadings.  In  static 
tests,  tunnel  closures  are  monitored  with  increasing  load,  pro¬ 
viding  a  continuous  closure-load  relationship.  It  is  also  recog¬ 
nized  that  in  dynamic  tests,  tunnel  closure  is  measured  only  after 
the  load  is  removed.  This  procedure  provides  only  one  closure- load 
relation,  i.e.,  the  final  tunnel  closure  versus  peak  dynamic  load. 
Therefore,  the, dynamic  closure  measurement  docs  not  represent  the 
peak  closure  which  may  differ  from  final  closure  due  to  recovery 
that  occurs  during  unloading.  During  unloading  the  medium  near 
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the  tunnel  may  also  experience  inelastic  behavior  (loss  of  strength 
in  low  pressure  states) ,  causing  tunnel  closure  measurements  to 
deviate  substantially  from  the  peak  closure. 

A  series  of  numerical  analyses  were  performed  to  provide 
information  regarding  the  importance  of  the  following  effects: 

(1)  Recovery  during  unloading 

(2)  Inertia  effects 

(3)  Material  rate  effects 

The  results  of  this  numerical  investigation  are  summarized 
in  Figure  6.6.  The  solid  line  illustrates  the  differences  between 
peak  closures  and  closures  measured  after  unloading  {1%  peak 
versus  S°0  .  This  result  was  generated  by  a  static  analysis  using 
nominal  material  properties.  The  dashed  line  shows  the  influence 
of  inertia  on  closure.  (As  anticipated,  the  inertia  effects  are 
not  significant,  i.e.,  SRI  dynamic  tests  are  essentially  quasi¬ 
static.)  The  dash-dot  line  shows  the  effect  of  increasing  the 
yield  strength  to  account  for  strain-rate  dependent  material 
properties.  Finally,  the  dash-double-dot  line  gives  the  effect  of 
decreasing  the  yield  strength  to  account  for  pore-water  pressure 
under  static  loading. 

The  increase  in  strength  is  assumed  to  account  for  strain  rate 
effects  (a  25%  increase  in  failure  strength  for  pressure  levels 
at  12  ksi)  .  Including  strain  rate  effects  and  recovery  due  to 
unloading  accounts  for  a  significant  portion  of  the  differences 
that  occur  between  tests  (for  =  12  ksi:  dynamic  closure  =  3.5',', 
static  closure  =  7%).  These  results  agree  very  well  with  the 
experimental  data. 

In  summary,  it  appears  that  material  strain  rate  is  the 
dominant  effect  in  reducing  tunnel  closure  in  dynamic  tests  and 
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Figure  6.6.  Comparison  of  Inertia  and  Strain  Rate  Effects  for 
Dynamic  Hydrostatic  Loading  Simulations. 
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porewater  pressure  is  the  dominant  effect  in  increasing  tunnel 
closure  in  static  tests.  Clearly  the  recovery  during  unloading 
is  appreciable,  but  secondary.  Inertia  effects,  as  expected,  are 
small . 

Because  of  the  substantial  differences  that  occurred  in  static 
and  dynamic  tests,  it  would  be  prudent  to  first  perform  static 
and  dynamic  material  tests  on  solid  material  specimens  (i.e.,  no 
tunnel)  in  future  test  series.  These  tests  would  indicate  material 
rate  and  porewater  effects  as  well  as  provide  a  data  base  regard¬ 
ing  strain  rate  and  porewater  effects. 
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SECTION  7 


SENSITIVITY  OF  TUNNEL  RESPONSE  TO 
LATERAL  CONFINEMENT 

In  this  section  the  results  of  a  series  of  numerical  analyses 
are  used  to  evaluate  and  interpret  the  sensitivity  of  small-scale 
tests  (i.e.,  tunnel  closures  and  lateral  boundary  behavior)  to 
over-  and  underconfinement.  An  explicit  definition  for  neutral 
confinement  of  a  finite  test  specimen  is  difficult  to  establish, 
although  four  reasonable  definitions  are: 

(1)  Uniform  lateral  confinement  pressure  that  corresponds 
to  the  result  obtained  for  tunnel-free  uniaxial  strain 
conditions.  This  gives  a  lower  bound  to  the  confinement 
pressure . 

(2)  Uniaxial  strain  (roller)  lateral  boundary  which  could 
result  in  a  non-uniform  lateral  confinement  pressure. 

This  gives  an  upper  bound  to  the  confinement  pressure. 

(3)  Uniform  lateral  confinement  pressure  required  to  maintain 
zero  lateral  displacement  (hoop  strain)  on  the  specimen 
at  the  location  corresponding  to  the  SRI  simulated*  uni¬ 
axial  test  (i.e.,  ~2a  above  the  mid-height). 

(4)  Uniform  lateral  confinement  pressure  required  to  maintain 
zero  lateral  displacement  (hoop  strain)  on  the  specimen 
at  the  location  corresponding  to  free-field  analyses 
(e.g.,  ~5a  above  the  mid-height  as  indicated  in  Section  5). 

Of  course,  as  the  specimen  boundary  increases  (relative  to  the  tun¬ 
nel  diameter,  2a)  all  four  definitions  become  more  equivalent.  It  is 
the  possible  interaction  of  the  tunnel  response  with  the  lateral 

*In  the  SRI  simulation  the  sum  of  two  strain  gauge  (180°  apart) 
measurements  is  set  to  zero.  This  usually  implies  that  one 
measurement  is  positive  and  the  other  is  negative,  and  not  that 
the  lateral  displacement  is  actually  zero. 
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boundary  behavior  which  makes  an  explicit  definition  for  neutral 
confinement  difficult. 

A  series  of  five  separate  analyses  were  performed  to  evaluate 
the  effects  of  uniform  lateral  confinement  pressure  on  tunnel 
response.  Figure  7.1  indicates  the  various  ratios  of  vertical 
confinement-lateral  confinement  pressures  usrjd  in  this  study.  The 
lateral  confinement  pressure  corresponding  to  a  tunnel-free  uniaxial 
strain  test  was  used  as  the  lower  bound  case.  The  lateral  confine¬ 
ment  pressure  was  increased  in  ten  percent  increments  for  the  other 
four  analyses.  The  tunnel  crown  and  springline  closures  obtained 
for  these  analyses  are  shown  in  Figure  7.2  and  7.3,  respectively. 
These  results  clearly  indicate  the  sensitivity  of  tunnel  crown 
closure  to  small  changes  in  lateral  pressure,  especially  near  lower 
bound  (underconfined)  conditions. 

The  associated  sensitivity  of  the  response  of  the  specimen's 
lateral  boundary  is  shown  in  Figure  7.4.  Note  that  as  the  confine¬ 
ment  pressure  is  increased,  the  lateral  bulging  decreases  with  zero 

* 

lateral  displacement  at  z  =  5a  for  Pu  =  -1.1P  and  at  z  =  2a  for 
P^  =  ~1.2P^  (where  P^  is  the  lower  bound  tunnel-free  confinement 
pressure) . 

In  order  to  show  the  sensitivity  of  tunnel  closure  to  confine¬ 
ment,  the  results  of  the  present  study  are  re-plotted  in  Figure  7.5 
and  compared  with  the  four  definitions  of  neutral  confinement. 
Examining  Figure  7.5  shows  tunnel  closures  to  be  extremely  sensitive 
to  underconfinement  using  all  four  definitions  of  neutral  confine¬ 
ment.  These  results  also  show  that  the  definition  of  neutral 
confinement  can  significantly  influence  the  behavior  of  tunnel 
closure  (and  lateral  boundary)  response.  On  the  other  hand,  as  the 
test  specimen  is  overconfined,  the  tunnel  crown  closure  approaches 
that  for  the  ideal  uniaxial  strain  (roller)  boundary  state.  Note, 
however,  that  the  overall  tunnel  closure  (e.g.,  springline)  response 
is  more  sensitive  to  overconfinement.  The  isotropic  (hydrostatic) 
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Figure  7.1.  Ratios  of  Applied  Stress/Confinement  Stress  Used  for  Sensitivity  Study 


VERTICAL  CLOSURE  (PER  CENT) 


(PER  CENT) 


cal  Boundary 


Lateral  Confinement  Pressure  Sensitivity  Effects  on  the  Strain  Profile 
at  the  6.4a  Boundary  for  a  Vertical  Applied  Pressure  of  8  ksi. 


EFFECTS  8F  LATERAL  CONFINEMENT  STRESS  ON  CRONN  CLOSURES 
VOF/SHXP 


loading  condition  gives  smaller  (less  conservative]  tunnel  closures 
than  uniaxial  strain  loading. 

The  numerical  calculations  regarding  the  sensitivity  to  lateral 
confinement  of  tunnel  closure  and  lateral  strain  at  the  mid-height 
(i.e.,  the  degree  of  lateral  bulging)  are  compared  with  the  SRI 
experimental  results  in  Figures  7.6  and  7.7,  respectively.  The 
experimental  closures  are  presented  for  lateral  confinement  pres¬ 
sures  corresponding  to  the  SRI  simulated  uniaxial  strain  condition 
(Definition  3  above) ,  20  percent  overconfinement  and  20  percent 
underconfinement.  The  theoretical  calculations  agree  very  well  with 
experimental  data  with  regard  to  the  sensitivity  of  closure  to 
lateral  confinement,  although  numerical  closure  results  are  con¬ 
sistently  smaller. 

Figure  7.7  compares  experimental  and  theoretical  lateral 
boundary  strains  at  mid-height.  Again,  the  theoretical  and  experi¬ 
mental  results  indicate  the  same  trends,  although  the  analytical 
model  predicts  lateral  motions  greater  than  the  experiment,  e.g., 
for  uniaxial  loading  at  =  7.5  ksi,  the  experimental  mid-height 
strain  was  .02$  compared  to  the  analytical  results  of  .08$,  a  factor 
of  4  greater. 

It  should  be  noted  that  the  SRI  data  showed  considcrab le 
scatter  in  lateral  strain  measurements  apparently  for  identical 
tests.  For  example,  six  simulated  uniaxial  strain  tests  measured 
a  scatter  in  lateral  strain  at  mid-height  for  a  vertical  pressure 
of  7.5  ksi  from  +.04  percent  (bulging)  to  -.01  percent  (inward 
boundary  motion).  Although  the  calculated  results  do  not  lie  with¬ 
in  ti\e  band  of  the  experimental,  the  numerical  calculations  do 
reflect  the  actual  bulging  observed  in  the  experiment.  Because  of 
the  scatter  observed  in  experiments  and  the  relative  sensitivity 
of  the  lateral  bulging  to  material  variations,  these  differences 
in  experimental  and  theoretical  results  should  not  be  unexpected. 
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LATERAL  STRAIN  (PERCENT) 


LATERAL  STRAIN  V  SAX  ANA  2C2 

COMPARISON  OP  TEST  OATA  AMO  ANALYTICAL  RESULTS 


Figure  7.7, 


VERTICAL  PRESSURE  (KSI) 


Comparison  between  Experimental  and  Theoretical 
Specimen  Lateral  Strain  at  Mid-Heigbt  Versus 
Vertical  Pressure,  for  Uniaxial  Strain,  Over¬ 
confined  and  Underconfined  Loadings  of  SRI  RMG 
2C2.  (Aluminum  liner,  a/h  =  11.5.  P^x  is 
lateral  pressure  required  to  maintain  uniaxial 
strain  as  defined  by  SRI.) 


SECTION  8 


NUMERICAL  SIMULATION  OF  STATIC 
UNIAXIAL  STRAIN  EXPERIMENTS 

The  numerical  simulation  of  the  SRI  test  was  performed  using 
the  generalized  plane  strain  model  described  in  Section  3  and  the 
combined  variable  friction  angle  and  strain  hardening  Drucker- 
Prager  material  model  described  in  Appendix  A.  Comparing  the 
results  of  these  analyses  with  those  for  ideal  uniaxial  strain 
shows  the  following  discrepancies  and  similarities  that  occur  in 
the  two  tests  in  regard  to  tunnel  closures,  and  the  lateral 
boundary  conditions. 

The  influence  of  the  tunnel  on  lateral  confinement  pressure 
is  illustrated  in  Figure  8.1.  Shown  here  are  the  ratios  of  the 
applied-confinement  pressures  corresponding  to  the  free-field 
and  the  SRI  simulated  uniaxial  strain  test.  Based  on  these  results, 
the  tunnel  appears  to  increase  the  lateral  confinement  stress  on 
the  lateral  boundary  by  approximately  twenty  percent.  This  explains 
the  underconfined  analytical  results  obtained  in  the  preliminary 
analysis  (Section  4)  where  the  numerical  simulation  of  the  SRI 
uniaxial  strain  test  was  performed  by  essentially  applying  lateral 
pressures  corresponding  to  free-field. 

Tunnel  crown  and  springline  closures  obtained  from  the 
numerical  analysis  for  an  ideal  uniaxial  strain  and  simulated  uni¬ 
axial  strain  tests  are  shown  in  Figures  8.2  and  8.3,  respectively. 
These  numerical  results  indicate  that  for  practical  purposes,  the 
SRI  test  duplicates  tunnel  closures  corresponding  to  ideal  uniaxial 
strain  conditions.  Although  tunnel  closures  are  equivalent,  the 
deformation  behavior  on  the  lateral  boundary  differs  significantly 
(see  Figures  8.4  and  8.5).  The  analysis  of  the  SRI  test,  Figure 
8.5,  shows  significant  axial  variations  in  bulging  along  the 
lateral  boundary  caused  by  the  uniform  lateral  pressure  applied  to 
the  cylinder  boundary.  Note  that  the  lateral  motion  is  zeroed  one 
tunnel  diameter  above  the  mid-height  (similar  to  the  SRI  test). 
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Figure  8.1.  Effects  of  the  Tunnel  on  the  Confinement  Stress  at  the  6.4a  Boundary. 
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(P**  is  the  confinement  pressure  adjusted  until  the  lateral  strain 
is  zero  at  point  B  located  at  y  =  6.4a,  z  =  2. a.) 

Figure  8.3.  Comparison  of  Tunnel  Springline  Closures  for  Plane 
Strain  and  Generalized  Plane  Strain  Analyses  with 
Boundary  Conditions  Applied  at  A. 4a. 
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train  Distribution  at  the  6.4a  Boundary  for  the  Plane  Strain  Analysis  with 
5ller  Boundaries  at  6.4a. 


HORIZONTAL  9 TRAIN  (PERCENT) 

(The  confinement  pressure  P**  is  adjusted  until  the  lateral  strain  is  zero  at  point  B  located  at 
y  =  6.4a,  z  =  2a. ) 

Figure  8.5.  Strain  Distribution  at  the  6.4a  Boundary  for  the  Generalized  Plane  Strain  Analysis 
with  Stress  Boundaries  at  6.4a. 


The  degree  of  bulging  follows  tunnel  closure,  i,e.,  the  closure 
and  lateral  bulging  dependency  on  applied  loads  are  similar. 

The  ideal  uniaxial  strain  analysis  indicates  that  the  lateral 
stress  P„  is  non-uniform ,  see  Figure  8 .  (> ,  an  effect  which  increases 
as  tunnel  closure  increases.  For  low  pressures,  (elastic  behavior) 
the  ideal  uniaxial  strain  results  produced  nearly  uniform  lateral 
pressures.  Despite  the  differences  in  the  lateral  stress  distribu¬ 
tions,  the  intensity  (average  lateral  stress)  of  the  lateral  forces 
are  nearly  equivalent  (compare  Figure  8.6  with  8.7).  This  sug¬ 
gests  that  tunnel  response  is  more  sensitive  to  the  intensity  of 
the  lateral  pressure  than  the  distribution . 

These  analyses  support  the  general  validity  of  the  SRI  experi¬ 
ments  as  a  means  for  predicting  tunnel  response  of  an  idealized 
uniaxial  strain  test  under  a  static  loading  environment.  However, 
some  of  the  results,  specifically  lateral  boundary  behavior  appear 
inconsistent.  For  example,  numerical  analyses  of  the  SRI  simulated 
uniaxial  strain  test  indicate  that  significant  lateral  bulging 
occurs  especially  at  high  pressure  levels.  This  leads  to  the 
obvious  question  as  to  how  the  SRI  test  appears  to  deviate  from  the 
ideal  uniaxial  strain  configuration,  yet  give  tunnel  closures  which 
are  essentially  equal.  The  numerical  investigation  (see  Section  7) 
on  the  sensitivity  of  tunnel  closure  and  boundary  behavior  to 
over-  and  underconfinement  suggests  that  lateral  boundary  behavior 
has  only  a  secondary  effect  on  tunnel  closures  when  overconfined.  Thus, 
any  simulated  uniaxial  strain  test  that  is  slightly  ovcrconfined 
will  essentially  duplicate  tunnel  response  of  ideal  uniaxial  strain 
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Figure  8.6.  Stress  Distribution  at  the  6.4a  Boundary  for  the  Plane 
Strain  Analysis  with  Roller  Boundaries  at  6.4a. 


APPENDIX  A 


COMBINED  VARIABLE  FRICTION  ANGLE  AND  STRAIN 
HARDENING  DRUCKER- PRAGER  MATERIAL  MODEL 

This  appendix  describes  the  fundamental  relations  for  elastic- 
plastic  materials.  The  mathematical  development  is  presented  for 
the  general  elastic-plastic  material  with  the  specific  equations 
given  for  the  combined  variable  friction  angle  and  strain  hardening 
Drucker-Prager  model.  Finally,  the  modeling  capabilities  of  the 
variable  Drucker-Prager  model  are  compared  with  material  test  data 
for  SRI  RMG  2C2. 

Incremental  elastic-plastic  behavior  (i.e.,  derivation  of 
the  incremental  stress  -  strain  relations)  requires  i  7 1  three  separate 
entities:  (1)  a  yield  criterion  that  specifies  three  dimensional 

states  of  stress  corresponding  to  initiation  of  plastic  flow; 

(2)  a  flow  rule  relating  the  plastic  strain  increments  to  the  cur¬ 
rent  stress  state;  and  (3)  a  hardening  or  softening  rule  which 
describes  how  the  yield  criterion  is  modified  with  the  accumulation 
of  plastic  strains. 

The  initial  yield  and  subsequent  yield  criterion  for  an  iso¬ 
tropic  hardening  formulation  can  be  written  as 

F (o.j,  WP)  =  0  (A- 1) 

where  wP  is  the  accumulated  plastic  work  defined  as 

eP. 

wP  =  /  a .  .de?  .  (A- 2) 

ij  ij 
o  J  J 

For  the  combined  variable  friction  angle  and  strain  hardening 
Drucker-Prager  yield  criterion  is  defin'd  as 
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and  a  and  k  are  material  parameters  which  are  functions  of  both 
mean  stress,  Ji,  and  plastic  work,  Wp 


a  =  a ( J  j ,  Wp) 
K  =  K ( J i ,  Wp^ 


(A-5) 


The  variation  in  a  and  K  with  respect  to  hydrostatic  pressure  and 
plastic  work  are  curve  fitted  to  match  experimental  data  using 
spline  functions,  <J> ,  where 


a  =  $  i  ( J  i )  4>2(WP) 
K  =  $  3  ( J  3  )  $40VP) 


(A-  6) 


The  plastic  strain  increments  were  related  to  current  stress 
state  through  the  use  of  an  associated  flow  rule.  With  an  associated 
flow  rule,  the  plastic  strain  increments  are  related  to  the  yield 
function  using  a  normality  rule,  i.e., 


( A  -  7 ) 
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For  the  variable  Drucker-Prager  model 
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Associate  flow  rules  are  particularly  advantageous  in  implicit 
codes  (e.g.,  NONSAP),  since  their  use  leads  to  symmetric  system 
stiffness  matrices.  Moreover,  in  static  analysis,  associated  flow 
rules  for  strain  hardening  materials  guarantee  uniqueness,  whereas 
non-associated  flow  rules  can  lead  to  non-unique  static  solutions. 

Having  established  the  yield  surface,  a  flow  rule  and  the 
hardening  characteristics,  the  incremental  elastic-plastic  stress- 
strain  relations  are  obtained  through  algebraic  manipulation.  From 
the  definition  of  the  yield  condition  (i.e.,  during  plastic  flow, 

F  -  0)  we  obtain 

dF  =  0  =  do  ♦  dWp  ( A-  9 ) 

a°ij  1J  3WP 


for  the  variable  Drucker-Prager  yield  condition 
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The  increment  plastic  work  in  terms  of  the  stress  and  the  plastic 
strain  increments  is  equal  to 


dWp  =  o .  . de?  . 

ij  ij 


CA- 11) 
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The  stress  increments  are  evaluated  from  the  matrix  of  elastic 
constants  and  the  elastic  strain  increments  given  by 
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Solving  for  X 
increments,  de 

stress  state) 
leads  to 


in  terms  of  the  current  stresses,  a--,  and  strain 


ij  ’ 

is  defined  in  terms  of  the  current 
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by  substituting  equations  A-7,  A-ll  and  A-12  into  A-9 
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(A- 13) 


The  incremental  stress- strain  relationships  are  finally  determined 
by  substituting  equations  A-13  and  A-7  into  A-12  which  gives  the 
following  tangent  constitutive  matrix 
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where 

fEl  3F  3F  El 
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3  F  pfil  3  F  3  F  3F 
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A  comparison  of  the  combined  variable  friction  angle  and  strain 
hardening  Drucker-Prager  model  with  the  actual  experimental  data 
provided  by  WES  is  shown  in  Figures  A- 1  through  A-3.  These  figures 
clearly  demonstrate  the  improved  modeling  capabilities  of  the 
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Figure  A- 1  .  Comparison  of  Failure  and  Yield  Surfaces  for  the  Combined  Variable  Friction  Angle  and  Strain 
Hardening  Versus  Per f ectly-Plastic  Drucker-Prager  Models  of  SRI  RMG  2C2.  (v=0.182,  <$>=2. 5°) 


PBtPECTLr  PUWTIC 


J _ I _ L _ I _ 1, 

0.2  0.3  0.4  0.5  0.6 

AXIAL  STRAIN  (PERCENT) 


Comparison  of  Unconfined  Compression  Stress-Strain 
Characteristics  for  the  Combined  Variable  Friction 
Angle  and  Strain  Hardening  Versus  Per fectly-Plastic 
Drucker-Prager  Models  of  SRI  RMG  2C2.  (E  =  1.16  x 
106  psl,  Ou  =  3200  psi.) 


ML  SIMM  APPLIED  STOSS  CONFINEMENT  STOSS 
HIP  VS.  retPECTLT  PUOTIC  OEUdEX  -  PEACES 


Figure  A-3.  Comparison  of  Uniaxial  Strain  Loading  Applied  Stress-Confinement  Stress  for 
the  Combined  Variable  Friction  Angle  and  Strain  Hardening  Versus  Perfectly- 
Plastic  Drucker-Prager  Models  of  SRI  RMG  2C2  (v  =  0.182,  (p  =  2.5°). 


combined  variable  friction  angle  and  strain  hardening  Drucker- 
Prager  model  over  the  perfectly-plastic  Drucker - Prager  (or 
Mohr-Coulomb)  model. 
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APPENDIX  B 


ELASTIC-PLASTIC  MODIFIED  MOHR-COULOMB 
MATERIAL  MODEL 

The  Drucker- Prager  criterion  originally  evolved  as  a  3-D 
generalization  to  the  planar  Mohr-Coulomb  yield  law.  The  Drucker- 
Prager  criterion  simplifies  the  mathematics  encountered  with 
associated  flow  rules.  In  fact,  the  3-D  Mohr-Coulomb  criterion 
cannot  be  directly  used  as  a  flow  rule  for  states  of  stress  having 
two  identical  principal  stresses  (e.g.,  triaxial  compression).  The 
Mohr-Coulomb  criterion  is  modified  to  remove  this  problem,  as 
detailed  in  the  following. 

The  Mohr-Coulomb  material  model  is  based  on  the  Coulomb  slip 
concept.  That  is,  the  shear  stress  (x)  required  for  simple  slip 
is  considered  to  depend  upon  the  cohesion  (c)  and  linearly  upon 
the  normal  stress  (a  )  to  the  slip  surface,  Figure  B-l.  This  can 
be  stated  as: 

x  =  c  -  o^tancj) 

where  <p  is  the  so-called  friction  angle.  In  terms  of  principal 
stresses  (oi>o2>a3)  this  becomes: 

oi  -  a3  =  2c  cos<p  -  (ai+o3)sin4> 

This  function,  which  can  be  considered  a  pressure  dependent  Tresca 
condition,  represents  a  hexagonal  pyramid  in  stress  space,  Figure 
B-2.  The  yield  function  depends  only  on  the  extreme  values  of 
stress  and  is  independent  of  the  intermediate  stress.  This  is 
shown  in  the  7r-plane  intersection  as  each  side  of  the  Mohr-Coulomb 
surface  is  parallel  to  a  stress  axis. 


Figure  B-2.  Mohr-Coulomb  Yield  Function. 
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To  put  the  yield  function  into  invarient  form,  it  is  conven¬ 
ient  to  use  a  set  of  stress  invariants  introduced  by  Nayak  and 

Zienkiewicz  111].  These  new  invarients  (a  ,  a,  6  )  are  related 

m’  ’  o 

to  the  conventional  invarients  (Jj,  J2,  J^)  by  the  following: 


-  i  Ji 


1  sin'1 
3 


.-r-  /3")  (  J  , /o  ) 


<  en  <  4 

6  —  0  —  6 


If  one  considers  the  n-plane  (i.e.,  the  stress  space  plane  per¬ 
pendicular  to  the  line  o1=o2=o3  and  defined  by  the  equation 
01+02+03  =  0)  located  at  the  origin  in  stress  space  (Figure  B-3), 
then  0q  serves  as  an  angular  measurement  from  the  line  of  pure 
shear  over  to  a  stress  axis  projection.  The  term  0  acts  as  a 
radial  measurement  outward  from  the  n-plane  origin,  so  0  and 
6q  act  as  "polar"  coordinates  in  the  Ti-plane. 

These  invarients  lend  themselves  well  to  the  Mohr-Coulomb 
yield  criterion,  which  can  be  expressed  as  a  yield  function,  F 
in  the  following  form: 

F  =  ocos0o  -  ( 1/ /3) os in0Qs incj)  -  c  cosi>  +  omsin<f> 

For  associated  flow  rules  (as  found  in  NONSAP)  plastic  flow 
occurs  normal  to  the  yield  surface,  i.e., 


deP.  = 

1J  30U 

where  deP.  and  A  are  the  incremental  plastic  strains  and  a  scalar, 

ij 


3 


normal 


(01=02=03) 

✓ 


Tr-Plane  (ai+a2+a3  =  0) 

03 


’"-l 


Figure  B-3.  The  fl-Plane  in  Stress  Space 


respectively.  Once  the  gradient  of  the  yield  surface  is  determined, 
the  incremental  stress  -  strain  relations  are  defined  using  the  pro¬ 
cedures  defined  in  Appendix  A. 

Using  the  Nayak  and  Zienkiewicz  invarients,  the  gradient 
becomes : 


+ 


3F 


30 


o) 


30 
_ o 

3d  J1, 
k 


or  simply  as: 


{q}  =  Cj  (q^)  +  C2  f q 2 d  +  {q^} 


whe  re 


{qi>  =  |  {1,1,1,0,0,0}T 

{q0}  =  {S  ,S  ,S  ,2x  ,2t  ,2t  }T 

M2  Za  x’  y’  z  yz’  xz’  xy 


where  S  ,  S  and  S  are  the  stress  deviators.  The  C  constants  for 
x '  y  z 

Mohr-Coulomb  are  given  as 
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sin<J) 


C 

C 

c 


1 

2 

3 


cos0  1  +  tanS  tan30  + 
o  o  o 

(/3  sin0Q  +  cos0Q  s int}> ) 

(2o2cos30  ) 
v  o 


sin^ 

/3 


(tan30o 


tan0  ) 
o 


A  problem  arises  when  the  Mohr-Coulomb  model  is  adapted  to  the 
numerical  modeling.  When  the  stress  field  is  such  that  any  two 
principal  stresses  are  equal,  the  stresses  lie  on  a  corner  of  the 
Mohr-Coulomb  yield  surface  (i.e.,  0q  =  t  ^-)  .  At  this  point,  there 
is  a  non-unique  normal  to  the  surface  and  the  gradient  calculations 
are  not  defined.  Numerically  these  calculations  "blowup"  (C^  -*•  » 
as  0Q  -  ±  J). 

To  correct  this  problem,  a  Modified  Mohr-Coulomb  yield  surface 
has  been  defined  which  is  identical  to  the  ordinary  Mohr-Coulomb 
law  except  near  the  vertices  of  the  polygonal  surface.  Here  the 
corners  have  been  rounded  off  to  a  circular  arc  within  a  certain 
tolerance,  ac,  of  the  corner  condition,  0Q  =  ±  ^  (Figure  B-4). 

The  gradient  vector  is  thus  fully  defined  everywhere,  including 
the  region  where  0Q  =  ±  Also,  the  gradient  of  a  stress  state 
on  the  yield  surface  is  now  a  continuous  function  as  the  stress 
state  traverses  from  one  flat  portion  of  the  surface  to  another. 

The  Modified  Mohr-Coulomb  yield  function,  F,  and  gradient 
constants  C^,  C^  and  C^  remain  the  same  as  Mohr-Coulomb  except 
on  the  rounded  corners,  where  they  can  be  expressed  as 

F  =  2o2  -  2/1  x  o  cosa  +  2/2  a  cosa  — -  tana 

I  2  c 

-  2  x  ip  tana  +  x2  ( 1  +  tan2a  ) 
r  2  C  r?;  C 


i  oy 


Ill 


_ t  /6~ _ 

o3(3  -  4sin2a) 

and  ac  is  the  angle  measuring  the  tolerance  region  for  rounding  at 

the  corners  (see  Figure  B-4)  .  (Note:  The  tolerance  angle  for 

0  =  need  not  be  the  same  as  that  for  0  = 

o  6  o  6  ' 

For  very  small  rounding  of  the  corners,  the  Modified  Mohr- 
Coulomb  model  gives  virtually  identical  results  to  the  ordinary 
Mohr-Coulomb  model  when  the  latter  model  can  be  used.  The  rounded 
surface,  however,  is  generally  applicable,  while  the  ordinary  Mohr- 
Coulomb  condition  is  not.  When  the  rounding  region  is  made  larger, 
the  Modified  Mohr-Coulomb  model  attains  a  form  similar  to  the  conical 
surface  in  the  Drucker-Prager  model. 

Figures  B-5  and  B-6  compare  the  differences  in  tunnel  closures 
that  occurred  with  Drucker-Prager  and  Modified  Mohr-Coulomb 
material  models.  In  all  cases,  the  Modified  Mohr-Coulomb  gave  more 
closure  than  the  Drucker-Prager.  This  was  expected,  since  the 
Drucker-Prager  material  constants  were  selected  so  that  the  Drucker- 
Prager  surface  conservatively  bounds  the  Modified  Mohr-Coulomb 
surface.  Nevertheless,  both  models  predicted  similar  trends  and 
differences  in  closure. 
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SRI  RMG  2C2  Closure  Comparison  Between  Drucker-Prager  and  Modified 
Mohr-Coulomb  for  Hydrostatic  Loading  with  Both  Roller  and  Stress 
Confinement. 


bRI  RnG  zCz  Closure  Comparison  Between  Drucker-Prager  and  Modifie< 
Mohr-Coulomb  for  Plane  Strain  Roller  Boundary  Conditions. 
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